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ABSTRACT 

We study the evolution of the ultra-relativistic shock wave in a plane-parallel atmosphere adjacent 
to a vacuum and the subsequent breakout phenomenon. When the density distribution has a power 
law with the distance from the surface, there is a self-similar motion of the fluid before and after the 
shock emergence. The time evolution of the Lorentz factor of the shock front is assumed to follow a 
power law when the time is measured from the moment at which the shock front reaches the surface. 

The power index is found to be determined by the condition for the flow to extend through a critical 
point. The energy spectrum of the ejected matter as a result of the shock breakout is derived and 
its dependence on the strength of the explosion is also deduced. The results are compared with the 
self-similar solution for the same problem with non-relativistic treatment. 

Subject headings: gamma rays: bursts — hydrodynamics — relativity — shock waves — supernovae: 
general 


1. INTRODUCTION 

In a supernova (SN) explosion of a compact star, 
a shock wave is accelerated to a relativistic speed 
in the surface layer with a sharp density gradi- 
ent and results in relativi s tic m o tion of the ejecta 
llG olgatc. McKee, fc Blevins! 1973 iMatzner fc McK'^ 
Il99'9t iTa.n. Matzner Kr. McKef This kind of SNe 

is observationally classified into type Ic (SNe Ic). From 
radio observat i ons f or one of SNe Ic (SN 1998bw), 
iKiilka.rni et inferred that the radio shock front 

originated from this SN was moving at relativistic speeds 
with the Lorentz factor between 1.6 and 2. Since 
the discovery of a close associati on of a T-ray burst 
GRB 980425 with t his supernova llGala.ma, et 31.1179981: 
llwamoto et al.l[T9^ . there has been accumulated evi¬ 
dence that some of SNe Ic a r e rela ted to a class of 7- 
ray bursts. iHees fc MeszarosI i|1992ll proposed a plausi¬ 
ble model (the fireball shock model) to describe 7-ray 
bursts as phenomena taking place at cosmological dis¬ 
tances, in which relativistic motion with the Lorentz fac¬ 
tor of ~ 100 is required to accommodate the observed 
variability on very short time scales to huge 7-ray in¬ 
tensities and the shock converts its kinetic energy of the 
baryons into electrons to emit nonthermal 7-ray photons. 
Other models involving no relativistic motion would lead 
to a conflict that the size of the source indicated by the 
observed short time scales are compact enough for 7- 
ray photons to collide to yield e+e“ pairs and emit little 
7-ravs. This model pre dicted afterglows of 7-ray bursts 
(lM4szaros fc Hee.j1 99711 ■ The first detection of the X-ray 
afterglow from GR B 970228 by the Beppo-Sax satellite 
ijGosta, et al.ll199^ and the subsequent identification of 
the host galaxy confirme d the 7-rav burst as an e vent at 
a cosmological distance llvan Paradiis et aT ][l993). Thus 
it will be important to investigate under what conditions 
the matter ejected from a SN Ic (or a core collapse SN 


in general) give rise to this relativistic motion. 

Another aspect of the relativistic SN ejecta can 
be seen in the production of light elements. Re¬ 
cent observations for metal-poor stars have shown 
that the abundance of Be increases proportionally to 
the abundance of Fe, which indicates that primary 
processe s predominantly produce Be. It has been 
shown llFields. Daigne, C^^^e. fc Va,ngioni-F1a,ml I2flfl2l 
lNa.ka,miira Shigeva,mall2flfl4ll that the relativistic mo¬ 
tion of ejecta from SNe Ic may contribute to the pro¬ 
duction of light elements Li, Be, and B through spal¬ 
lation reactions as a primary process because the sur¬ 
face layer of the ejecta is composed of G and O. 
Though these studies are based on spherically symmet¬ 
ric models for SNe Ic, it would be interesting to in¬ 
vestigate the effects of an aspherical SN explosion on 
the acceleration of the ejecta and the productio n of 
the light elemen ts. iTan. Matzne r fc McKS 3 ll2001fl and 
iNakamura fc Shigevamal ll2004ll pointed out that an as¬ 
pherical explosion may enhance the acceleration com¬ 
pared to a spherical explosion with the same explo¬ 
sion energy. Thus an aspherical explosion might pro¬ 
duce more light elements. Furthermore, the acceleration 
strengthened by an aspherical explosion should be real¬ 
ized by a collimation of the flow. The collimated flow or 
jet is required to interpret th e time evoluti on of the af¬ 
terglow of some 7-ray bursts ||Rh oa dimi^). Hence the 
connection between SNe and 7-ray bursts needs to be 
explored in this respect as well. A linear stability anal¬ 
ysis of the corresponding spherical or plane-parallel flow 
will be a good starting point to understand this possible 
collimation mechanism. 

To study the acceleration of matter before and af¬ 
ter the passage of a shock wave across the surface of 
a stratified medium adjacent to a vacuum in general, we 
will present a self-similar solution for an ultra-relativistic 
shock wave propagating in a plane-parallel medium. This 
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problem with out taking into account relativistic effects 
was posed by ICandel’m an fc Fra,uk-Ka ,menetskiil ill 95fif) 
and fully solved bv ISakurail (119601) . Their solutions in¬ 
evitably involve an infinite velocity at the surface after 
the shock emergence. Therefore a self-consistent solution 
needs to take into account relativistic effects. To obtain a 
relativ istic solution, we follow the procedure of iSakur^ 
KTml) with a different choice of n on-dimensional inde¬ 
pende nt variables. As was done in IBlandford fc McKeel 
(Il97(ill for an ultra-relativistic shock wave emanating 
from a strong point explosion, the Lorentz factor of the 
shock front is assumed to evolve with time For 

our problem, the time t is assumed to be measured from 
the moment when the shock front hits the surface and 
the parameter m is to be determined by the condition 
that the flow extend through a critical point (see Q . 
A similar work was studied bv iPerua Kr. Vietril ll2nn2ti to 
describe the propagation of a relativistic shock wave in 
an exponential atmosphere. In this situation, the flow 
after the shock emergence cannot be treated because the 
surface is located at infinity. Another work related t o 
the present work was done bv l.Tohnson fc McKeel lll971ll . 
They derived a solution for free expansion into a vacuum 
with the initial condition that the fluid is at rest while we 
will be concerned with a solution for f ree expansion of the 
shocke d fluid moving into a vacuum. I.Tohnson fc McK^ 
liTMIl also derived a solution for the shock propagation 
into a medium of decreasing density which starts from 
two adjacent fluids at rest with different pressures. They 
have treated the same problem as in the present paper 
except for the initial conditions. Thus the solutions will 
be compared in some detail. 

Since the assumed density distribution is realized when 
a polytropic gas is in hydrostatic equilibrium with a con¬ 
stant gravitational acceleration, there will be many situ¬ 
ations to which this solution can apply in addition to the 
surface of a star. The analytical description of the flow 
presented in this paper will enable us to easily investi¬ 
gate the linear stability of the flow. Therefore this study 
is essential to investigating whether a disturbed shock 
front will enhance the acceleration of matter and lead 
to a collimated flow. Furthermore, this study will give 
a possible physical mechanism of the energy injection to 
the fireball model. 

Next section outlines the problem discussed in this pa¬ 
per. The flow before shock emergence is formulated and 
its self-similar solution is derived in m The subsequent 
self-similar flow after shock emergence is described in m 
The analytical solution for the flow at extremely large t 
is also presented. Based on these solutions, discusses 
the energy spectrum of the outermost layer of the ejected 
matter as a result of the shock breakout into a vacuum. 

concludes the paper. The Appendix gives the ana¬ 
lytical description for the self-similar flow before shock 
emergence. 


2. OUTLINE OF THE PROBLEM 

A self-similar evolution of the flow associated with the 
ultra-relativistic shock wave in a plane parallel stratified 
medium and the subsequent flow after the shock emer¬ 
gence is investigated in the following context: 

(a) The initial density distribution po{x) is given by a 


power law as 
Po{x) 


( bx^ for x > 0 , 
\o for X < 0, 


( 1 ) 


where x is the distance measured from the surface. 
Two positive constants b and <5 have been intro¬ 
duced. Note, however, that this distribution is a 
good approximation only near the surface of a star 
because the gravitational acceleration needs to be 
assumed constant and the sphericity is ignored. We 
will consider two specific cases in some detail: a 
convective envelope expressed with (5 = 1.5 and 
radiative one with (5 = 3, though our solutions pre¬ 
sented in this paper hold for any positive 6. The 
exterior of the star denoted by negative x is as¬ 
sumed to be vacuum before the shock emergence. 


(b) An ultra-relativistic shock wave propagates toward 

the surface in a static medium with the density 
distribution given above. The effect of gravity on 
the shock wave is assumed to be negligible. The 
pressure in front of the shock wave is neglected. 
The time t is measured from the instant of shock 
arrival at the surface. Thus this stage before the 
shock emergence is denoted by negative t. 

(c) After the shock emergence (t > 0), the fluid expands 

into the vacuum. The initial conditions for the flow 
are given by the solution in the stage (b) when the 
shock hits the surface {t —0). The analytical 
expression for the asymptotic solution with t —>■ oo 
will be shown. 


3. BEFORE SHOCK EMERGENCE 
3.1. Formulation 

In this section, we derive equations governing the prop¬ 
agation of an ultra-relativistic shock wave and the associ¬ 
ated self-similar flow from basic equations of relativistic 
hydrodynamics and discuss the boundary conditions at 
the shock front. Next, we define the similarity variable 
and convert the partial differential equations into a set 
of ordinary differential equations. The energy and mo¬ 
mentum conservation laws in a fixed frame (x, t) can be 
expressed as 

^(7^(e-P/3^p))-f ^(7^/3(e-Pp)) = 0, (2) 


f) f) ciT) 

^(7^/3(e+.)) + ^(7^/3^(eTp)) + £=0, (3) 

in a plane parallel atmosphere. Here 7 denotes the 
Lorentz factor, /3 the velocity divided by the speed of 
light c (we set c = 1 in the following), and p denotes the 
pressure. The energy density e per unit volume includes 
the thermal energy and the rest mass energy. The conti¬ 
nuity equation is written in terms of the number density 
n' in the fixed frame as 


dn' 


-(, 5 „')= 0 . 


( 4 ) 


This density is defined by n' = ny with the number den¬ 
sity n in the co-moving frame. Equations ©■^0 are 
to be solved with an equation of state. Since we are 
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concerned with the propagation of a shock wave in the 
relativistic limit, it will be appropriate to use the ultra- 
relativistic equation of state 

e = 3p. (5) 

Using the relation between the partial time derivative 
d/dt at a fixed point and the total derivative d/dt'. 
d/dt = djdt -\- pd/dx, together with this equation of 
state, equations @-0 are combined to yield 

;i = _4g, ,7, 

= ( 8 ) 

Then the non-dimensional pressure /, Lorentz factor 
g, and density h are introduced as functions of the sim¬ 
ilarity variable x- They are defined to become unity at 
the shock front satisfying th e ultra-relativistic Rankine- 
Hugoniot relations (see e.g., iBlandford fc McK^IlQTlH 
as follows, 

p = ^6XfrV(x), (9) 

7" = (10) 

n' = 2-xfr^hix), ( 11 ) 

g 

where Xg and P are the position and the Lorentz factor of 
the shock front and p is the mean mass of the constituent 
particles. As was mentioned in §1, the Lorentz factor of 
the shock front P is assumed to evolve with the time t as 


These relations hold to the accuracy of the order of 
P“^. Using these relations and substituting Pl-(fTHl into 
equations 10-0 yields the following ordinary differen¬ 
tial equations. 


Idlnf (S — m)gx + 8m — AS 
g dx {m + l) {g'^x^ - 8gx + 4) ’ 


(18) 


Idhig —mgxA-7m — 8d 

g dx im + l)ig‘^x^ - 8gx + A)' 


(19) 


1 din ft, (m — S)g^x'^ + ( 81 J — IOto)^^ -|- 18to — lOd 
g dx (m-bl)( 2 -gx)( 5 ^X^- 85 X + 4) 

( 20 ) 

These equations have been expanded in powers of P ^ 
and 7 “^ and truncated at the first contributions. Given 
values of 5 and to, these equations are integrated from 
X = 1 to X = —00 with the boundary conditions jnj. 
However, it is found that the derivatives diverge at a fi¬ 
nite X except for a certain value of to. To avoid this, the 
denominators and numerators on the right hand sides of 
these equations must vanish simultaneously, which oc¬ 
curs when 

TO = (t2V3-3)(5. (21) 

The shock being accelerated toward the surface indicates 
a positive to (= [2^/3 — 3) 6). If the negative value were 
chosen, the shock Lorentz factor would decrease to zero 
(less than 1) as the shock approaches the surface. This 
is not physical and the approximation of the relativistic 
limit would not hold. Hence, equation m is reduced to 
a simpler form, 

1 ding ^_V35_ 

g dx {2 + v^(J+l)}(gx-4-2y3)' ^ ^ 


( 12 ) 

where A is a constant. Following IHI and ford fc McKeel 
dia, the similarity variable x is defined as 

X = {l + 2(TO + l)Pn(l + f). (13) 

Here the power index to must be greater than —1 so 
that the shock speed not exceed the speed of light. Note 
that when the shock reaches the surface, i.e. t — 0 , 
the Lorentz factor P diverges, i.e., the shock speed ap¬ 
proaches to the speed of light. The value of to can not 
be determined from the initial conditions, but from the 
condition that the solution extend across a critical point 
to infinity. The boundary conditions at the shock front 
are expressed as 


/(l) = g(l) = ft(l) = l. (14) 

The coordinate system can be transformed from (a;, t) 
to (P, x) by the following relations 

i^ = 2(m + ir|-, (16) 

d toP d , _,,,2 . d 

'S = -—+ + 


The o ther t wo eq uations are also simplified to equations 
(fTT^ and iTHTt in the Appendix . It should be noted 
here that I.Tohnson fc McKea lll971ll have also presented 
an expression for the shock propagation in equation (43) 
of their paper, which indicates 

P oc n 2 ( 2 +V 3 ), (23) 

This relation is derived from an invariant quantity car¬ 
ried by forward characteristics and should hold in gen¬ 
eral situations. Since the density n at the shock front 
evolves with time t as n cx this expression is pre¬ 

cisely equivalent to the shock propagation indicated by 
the value of to in equation 

3 . 2 . Solutions 

Though the analytical form obtained by integrations of 
HBj-innii is available and presented in the Appendix, the 
dependent variables /, g, and ft can be expressed in terms 
of gx rather than the independent variable x- Thus we 
numerically integrate equations to obtain the 

distributions of the pressure, density, and Lorentz factor 
as functions of x- The results are shown in Figure 1 
for (5 = 3 and Figure 2 for 5 = 1.5. This expression 
is more convenient for investigating the distribution and 
time evolution of the flow. 

Here we will m ake a comparison of this s olution with 
that presented bv I.Tohnson fc McKeel lll971ll . They have 
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Fig. 1.— The solution of equations fR-m for (5 = 3 which 
corresponds to m = 6\/3 — 9 ~ 1.392 from equation J21J. The inte¬ 
gral curve g passes through the singular point which is represented 
by gx = 4 — 2\/3 without any special behavior. 



Fig. 2.— The same as Figure^but for S = 1.5 which corresponds 
to m = 3V3 - 4.5 ~ 0.696. 


(derived a relation between the Lorentz factors 7 and the 
positions a; of a fluid element at two different moments 
after passing through the ultra-relativistic shock front 
propagating in the medium with the density distribution 
given by equation o, 

xo 2{l-by3(l + (5)} 

(24) 

Here 70 and aio denote the values at the moment when 
the fluid elemen t is hit by the shock wave . Note that the 
original form in I.Tohnson k McK^ lll971ll is slightly dif¬ 
ferent from the above expression because our coordinate 
X is different from theirs. We have confirmed that this 
relation can be deduced from our solution by the follow¬ 
ing procedure. In the relativistic limit, the distance x 
from the surface is simply proportional to time t. Thus 
we obtain 


X t 

Xo to 




(V3+2) 

•/35 


(25) 


Here the time t has been converted into the shock Lorentz 
factor r using equation OS- Then the adiabatic evolu¬ 
tion of the fluid element after passing the shock front 
yields a relation between the non-dimensional functions 


and the Lorentz factors, 


gsvsfh 


-3 

4 


4(1 + V3) 
^ 70 ^ 3^3 


(26) 


Elimi n ating g, /, h from this equation and equations 
IIAI 8 II (IA20II in the Appendix, gx can be expressed as a 
function of the ratio of the Lorentz factors 70 / 7 . Substi¬ 
tuting this expression into ^ = —gx in equation ifTmli . 
then g is expressed as 


Ino- {3-k273(l-k(5)} A - 1 

V3+2 2{1 + V3(1 + ( 5 )}A ’ 

where A = j . (27) 

Eliminating g from equations (P|l and m yields the re¬ 
lation ll^ . Thus the assumption of self-similarity does 
not lose any generality. The above relation (EH) can yield 
the ratio of the Lorentz factors 70/7 of a fluid element 
when it passes the original surface by setting x = Q. 
At this moment, the shock is already emerging from the 
surface. On the other hand, the self-similar solution pre¬ 
sented in this section holds as long as the shock exists. 
To describe the flow after the shock disappears as it hits 
the surface, an extended solution is developed in the fol¬ 
lowing sections. 


3.3. Limiting Behavior 

In order to extend the solution obtained in the previous 
section to describe the flow after the shock emergence, 
the limiting behavior of the fluid when t —> —0 needs to 
be described. From equation EH, the position of the 
shock front is given by the formula 

2 (TO-M)r 2 }' 

Hence Xg ^ —t as t ^ —0. Thus the behavior of /, g, 
and h defined by (©-(EH is determined by the condition 
that p, 7 , and n' must not vanish or diverge at all places 
with finite x’s. To satisfy this condition, we need to 
eliminate their dependences on t to yield 

8-m 2(5 

/i(x) oc (-X)^ cx (-x) 2+'^(i+'5), (29) 

-m 

5 i(x) cx (-x)^ fx (-x) =+ 3 ^( 1 +^), (30) 


S — m 2(5 

^i(x) fx (-x)^ fx (-x)2+'^(i+3), (31) 


where the subscript 1 indicates the limiting functions. In 
the Appendix, this argument is verified in the analytical 
form. As a result, their distributions at the instant of 
the shock emergence are given by 


pi{x) = 


2Ab 


2A{l-h(2V3-3)5} 
3-k 2^3(1+ 5) 


2+V3(l + S) 


^ 2+\/3{l + (5) _ 


(32) 




A 


12 + 7^3-k 2 (2V3 + 3) (5 
2 [2-h v^(l + (5)] A 


_ \/3S _ 

2+V3(l + S) 


V3S 

X 2+V3(l + 5) ^ 


(33) 
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n[{x) 


2Ab 


2A{1 + ( 2 V 3 - 3),5} 

v^(l + ^) 

{3-h2V3(l-h(5)}Hi+v^(i+^)} 


25 

2+V3(1 + 5) 

X 


(34) 

Here the subscript 1 indicates distributions at t = 0. 
The distribution of the density has a clearly different 
power law index from that in the initial condition. This 
is a purely relativisti c effect, since it does not happen in 
non-relativistic case llSakurailll960ll . in which the strong 
shock always compresses the fluid by a constant factor. 


4. AFTER SHOCK EMERGENCE 
4.1. Formulation 


Since no new scale is introduced after the shock emer¬ 
gence, t he fluid keeps its self-similarity (in Newtonian 
case, see lSakura^l96^)^ . Lagrangian coordinates are con¬ 
venient for describing the motion of the expanding gas 
into a vacuum starting from given initial distributions 
llSakurailll96(ll . Thus we will use a Lagrangian coordi¬ 
nate a defined as it coincides with the Eulerian coordi¬ 
nate a; at t = 0 in the previous section. The relation 
between x and a is given by 


dx n\ 
da n' 

Now we seek a solution in the form of 


(35) 


p = Pi{a)F{vi), (36) 


7 ' = ll{a)G{p), (37) 


n = ni{a)H{r]), 


(38) 


X = aR{ri), 

where the similarity variable rj is defined as 


t 1 -c- 

rj = oc ta ^"+1 oc ta i+( 2 V 3 - 3)5 ^ 

7fa 

The initial conditions are 


(39) 

(40) 


E(0) = G(0) = iL(0) = R(0) = 1. (41) 


Using the following transformations from (t, a) to (rj, a) 
(here, we express the total derivative d/dt as the partial 
derivative d/dt in order to emphasize that it is taken 
with a =const.). 


d d 
^dt ^'^drj^ 


(42) 


d ij d ^ d 

da m + 1 dij da' 


(43) 


and 123), equations 13-® are converted to the following 
equations in the relativistic limit. 


G' 3F' H frjG' \ 

'G-~^~ (m + l)VG KG 

F' (5-m)H^/G-2(m+l)G' 
'F ~ (m + l)G + rjH^ ’ 


(44) 

(45) 


25 

2 + V3(l + d) ^ 



Fig. 3.— The solution of equations 1441- I4HI for <5 = 3. 



Fig. 4.— The same as Figure|21but for <5 = 1.5. 


H' _ 3F' 
ll ~ 4F ’ 

R' _ m + 1 / 1 \ 

'R~ rj Iv ~ RHVg) ’ 


(46) 

(47) 


where the prime denotes the derivative with respect to rj 
and m = (2^/3 — 3)<5 as in 123. 


4 . 2 . Solutions 

The ordinary differential equations 133-133 are nu¬ 
merically integrated with the initial conditions 1411) . Un¬ 
like the flow before shock emergence, the solution does 
not have a critical point. The results for 5 = 3 and 1.5 
are shown in Figures 3 and 4. 

There are two invariants in this solution. One is from 
the adiabatic condition (eq. I@ 3 )) 

Firj)H(rj)-F^ = 1, (48) 

and the other is 

G(7j)F(p)'^/^ = 1, (49) 

obtained by eliminating H{rj) from equations 14411 and 
133 - This is the Riemann inv ariant whose precise form i n 
the relativistic limit is given in i, Johnson fc McKeel 1197111 . 
From equations (EIi-iITtI). we can derive the asymptotic 
solution in the limit of 77 00 undergoing power law 

evolution : 

F{rj) = Girj) = G^rj^, H(rj) = H^rj^. (50) 
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Substituting them into equations 1011)-(EH), it follows 
that the power law indices must satisfy the following re¬ 
lations 


3 d , : N 

6 iF + m), 

2 TO -I- 1 

(51) 

(d — m)d — 2(to -I- l)ip 
m-\-l + d 

(52) 


(53) 

lb 

- - 6» = 1, 

2 

(54) 


where d = Hi/^/Gi. These equations have two solutions 
one of which is independent of (5 or m except for d. The 
other solution gives a negative ip which is unphysical. 
Thus the time scale of the asymptotic behavior does not 
depend upon the initial conditions because the leading 
edge of the fluid expands as fast as the speed of light. As 
a result, we obtain 

0 = 2 (^ - 1 ), = 73 - 1 , 0 = 


and d = 


2(2-73){737-H)-f2} 
1 -l- '\/3(<5 -l-1) 


(55) 


Using these power law indices, the non-dimensional po¬ 
sition R is found to evolve as from equation 07 |). 
We have confirmed that these are exactly matched with 
the results of the numerical integration. Furthermore, 
we can express the coefficients Fi , Gi and Hi in terms 
of S using equations TO and TO together with the ex¬ 
pression for d in iTOll as 

Fi = = 7-^, Hi = (56) 

These values also agree with our numerical results. All 
the coefficients depend on <5 and remember the initial 
density distribution. 

Our results are similar to those of l.Tohnson & McKed 
<ITm^ in some sense, though their initial conditions 
are different from ours. For example, it follows from 
our results that 'yj = 7i(Pi/P/)'^^^i where the sub¬ 
script / denotes the final value. On the other hand, 
I.Tohnson fc McKe'3 lll971ll show i'y'j = (pi/p/)^/^. 
Their fluid is initially static and uniform while ours is 
already not uniform and expanding at this initial mo¬ 
ment, though both treat a fluid expanding into a vacuum. 
Due to their initial conditions, their problem inevitably 
involves non-relativistic flow in which the Riemann in¬ 
variant obtained above in the relativistic limit is not ap¬ 
propriate. This causes the difference in the factor of 4 in 
front of 7 ^. Our approach can be applied to any small 
value of 6 as long as it is positive, because the shock front 
is eventually accelerated to ultra-relativistic speeds. It 
is, however, not ap propriate for 6 = 0 , a spe cial case of 
which is treated bv I.Tohnson fc McKed flTnl) . 


TABLE 1 

Power-law indices eor p and n . 


5 

t 

01np/51na 

relativistic non-relativistic 

a In 

relativistic 

n/d\n.a 

non-relativistic 

3 

0 

0.67 

1.89 

1.09 

3.00 


-l-CXD 

1.03 

3.96 

1.23 

4.56 

1 Fi 

0 

0.47 

0.84 

0.85 

1.50 


+ O0 

0.97 

2.54 

1.05 

2.79 


4.3. Asymptotic Distribution 

The distributions of physical quantities in the fluid in 
the limit ofrj—^oo are discussed in this section. The dis¬ 
tributions of our self-similar solution converge to power 
law distributions in this limit, which is referred to the 
asymptotic distribution here. This distribution is still 
intermediate and eventually changes to non-self-similar 
distribution when the matter cools down and no longer 
behaves as ultra-relativistic gas. Of course, a real star 
is spherical thus sphericity should be taken into account 
when we consider the motion far from the stellar sur¬ 
face. The effects of sphericity will be incorporated by 
expanding the flow variables in the power series of the 
ratio of the distance from the su rface to the stellar ra¬ 
dius llKazhdan fc Miirzinal 1 1993) or by introducing an 
empirical formula obtained from comparisons with nu- 
merical calculations for sp herically symmetric explosions 
llMatzner fc McK^ll999ll . 

When t oo, i.e. rj —>■ oo, the asymptotic dist ribu- 
tions of F, G and H have been already obtained in (Oil 
and (TO) . Thus the pressure, the Lorentz factor, and the 
density have power law distributions with respect to a as 


2(V35 + 1 + V3) 

P OC a 2y3+3(« + l) ^ 

(57) 

^ l+^3(5+l) 

7 OC a 2+y3(5+i) ^ 

(58) 

(4+vC3)5 + V3(l + V^) 

n cx a. 2{2+>/3(5+i)} ^ 

(59) 

Numerically, the pressure and the density have distribu¬ 
tions with respect to a as 

p OC n OC = 3), 

(60) 

p OC n OC = 1.5), 

(61) 


for specific values of 6. Comparing these results with 
those at the instant of shock emergence, there is only 
a small difference in the power law indices in the den¬ 
sity distributions as compared with the difference in the 
pressure distribution (see Table 1). 

In non-relativistic case, however, the distributions sig¬ 
nificantly change after the shock emergence. According 
to Saku rai’s solutions for a gas with the adiabatic in¬ 
dex 4/3 llSakurailll96nll . the distributions of pressure and 
density in the limit of t ^ oo are 

p oc n oc 7 ®®((5 = 3), (62) 

p oc n oc a^'™(i5 = 1.5). (63) 
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These power law indices are quite different from those 
at the instant of shock emergence (t = 0 ), details of the 
indices are shown in Table 1. 

The asymptotic behavior of G suggests that the 
Lorentz factor continues to increase without bound as 
long as the ultra-relativistic equation of state © holds. 
The adiabatic expansion will eventually invalidate this 
equation of state and cease a further acceleration. This 
effect will be taken into account to obtain the energy 
spectrum of the ejecta in the next section. 

5. ENERGY SPECTRUM 


From the self-similar solution discussed in the preced¬ 
ing sections, the energy spectrum in the ejected matter 
near the leading edge and its dependence on the initial 
strength of the explosion will be derived in this section. 
The initial strength of the explosion is deduced from the 
coefficient A in equation It is expected that A is 

proportional to the square of the explosion energy per 
unit mass, though both of the total energy and ejected 
mass involved in this self-similar solution are infinite. 

The energy per unit mass excluding the rest mass en¬ 
ergy is given by 



where 

T = {p + 4 p) 7 ^ — p — fry. (65) 


Here we have used the ultra-relativistic equation of state 

®. 

The energy spectrum is defined as the total mass with 
the energy per unit mass greater than e, 

5 + 1 

M{> e)= / fj,n[{x) dx (X {a{e)A)'^+^ ( 66 ) 

Jo 

(2+v/3)(5 + l) 

oc (a(e)i 4 ) 2 +v 3 {i+ 5 ) ^ 


where a(e) denotes the Lagrangian coordinate of the fluid 
ele ment with the energy per unit mass of e. Substitut¬ 
ing (for t = 0 ) or lf^ - (l^ (for t oo) into 

and taking the dominant term at a specific time, we 
can evaluate the relation between a and e. The energy 
distribution will be derived for limiting cases of t = 0 
and t = oo here and also compared with Sakurai’s non- 
relativistic cases for 6 = 1.5 and 3. 

At the instant of shock emergence, i.e., t = 0 , the 
specific energy e in (IMIl becomes dominated by the term 
p^/p near the edge of the ejected matter. Thus using the 
expressions for p, 7 , and p given by we obtain 


1 (m + l) 2+v^ {2 + v^(l + 5)} 

a(e) oc A^e oc A e 


(67) 


Then combining this result with equation yields 


M(> e) oc 



(5 + l)(2+V31 



( 68 ) 


For specific values of <5, the spectrum becomes 


M(> e) oc j 

r ^2.87^-2.87 

1 ^3.59^—3.59 

for 6 = 
for (5 = 

3, 

1.5. 

(69) 

On the other hand, Sakurai’s 
relativistic treatments gives 

solution from 

non- 

M(> e) cx j 

(A8/3e-4.37 

for S = 
for (5 = 

3, 

1.5, 

(70) 


where Anr represents the strength of the initial explosion, 
defined as 14 = AnrAj^, where 14 and Xs denotes the 
velocity and position of the shock front, A is t he eigen 
value determined by the method presented by iSakurail 
ifT^ . Thus this coefficient Anr is expected to scale 
with the square root of the explosion energy per unit 
mass. 

The distributions of F, G and H in the limit of t —> 00 
are given by pljl . g 21 ), and (P|l . In this limiting situa¬ 
tion, the internal energy has been completely converted 
to the kinetic energy. Thus 6 ^ 7 . Using these asymp¬ 
totic behavior, the same procedure leads to 


V3 —2 2(Tn + l) _ 1 _ _ 2{2 + y3(l + 5)]- 

a(e) cx oc Ai+Y3(i+'’) e i+%/3(i+5) ^ 

(71) 

and 


(2+v/3)(5 + l) 2(2 + V3)(5 + l) 

M(> e) OC A e . (72) 


Numerically, 


M(> e) cx 


r j^i.88^-3.77 

1 


for (5 = 3, 
for S = 1.5. 


(73) 


The smaller power index of e indicates that the ac¬ 
celeration is more effective in the inner part after the 
shock emergence. This phenomenon is not seen in non- 
relativistic treatments. In Sakurai’s non-relativistic so¬ 
lution, we obtain the energy spectrum with exactly the 
same dependence on Anr and e as at t = 0. This so¬ 
lution shows that the energy in the outermost layer of 
the ejected matter is dominated by the kinetic energy at 
both t = 0 and t ^ 00 and that the velocity profile in this 
region converges to one with the same power law index 
in both of these two limits. Therefore, it shows the fi¬ 
nal energy spectrum of the ejected matter w ith the same 
dependence on Anr and e as at t = 0 . Thus iFields et alJ 
ll2nf)l used the relation CZOl) for i5 = 3 to derive the em¬ 
pirical formula for the energy spectrum of the ejected 
matter M(> e) oc ® (Fex/ATej)^ ®, where the explo¬ 
sion energy per unit mass is expressed by the ratio of the 
explosion energy Fex to the ejected mass Mej. 

The above energy spectrum m as a result of the rela¬ 
tivistic shock breakout is still intermediate and the final 
energy spectrum should be modified from it to that at 
the moments when the acceleration ceases in each fluid 
element. Assuming that t he pressure suddenly dro ps to 
zero when e/p reaches 1 llJohnson fc McKeelH97I|l . the 
Lorentz factor 7 at this moment in each fluid element 
can be obtained as a function of a. Thus the final energy 
spectrum will change to 


(2+v^)(l + 5) (3+v^)(5 + i) 

M(> e) OC A v^'5 e 


(74) 


This leads to M{> e) oc e (Fex/ATej)^'^^ for 5 = 
3. For 5 = 1.5 the energy spectrum becomes M{> 
e) cx e“2.63 . On the other hand, the 

non-relativistic treatment yields the energy spectrum 
M{> e) cx {Ee^/Mgj)'^'* for 6 = 1.5 by substitut¬ 
ing Anr cx Y^Fex/Afej into the secon d equation of iTTnil 
which coincides with equation (24) of iMatzner fc McKed 
iTT^ . 
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6. CONCLUSIONS AND DISCUSSION 

We have presented a self-similar solution for a pla¬ 
nar flow associated with the ultra-relativistic shock wave 
propagating in a power-law density distribution adjacent 
to a vacuum. The solution can describe the flow before 
and after the shock front reaches the surface. The so¬ 
lution at the shock emergence is explicitly given in the 
analytical form as the initial condition for the flow af¬ 
ter the shock emergence. The asymptotic distributions 
of the flow as a result of the phenomenon are also pre¬ 
sented in the analytical form. 

Thanks to the det a iled c omparison with the results 
of I, Johnson fc McKeel lll971ll presented in 312 one can 
trace back each fluid element in the final state to the 
original position. Equation applied to the flow at 
t = 0 yields the ratio of the Lagrangian coordinate a to 
the original position xq as 

„ {3 + 2V3(1 + 5)}(^) 

xo~ 2{l-Lv^(l-f (5)} 

(75) 

These two coordinates are related to each other by the 
expression for 71 obtained in equation 03 in terms of 

the Lagrangian coordinate a and 7 g = /2 as 

indicated by equation m- Therefore the flow in the final 
state described in terms of a in 312 can be expressed as 
functions of Xq. 

It would be possible to formulate the flow before shock 
emergence with a Lagrangian coordinate as was shown by 


iMatzner fc McK^ lll999fl for non-relativistic flow. How¬ 
ever, it might not be possible to derive the eigen value 
of m analytically from this formalism because the resul¬ 
tant differential equation corresponding to equation ini 
would inevitably involve the self-similar variable for den¬ 
sity. Furthermore, it is not obvious how to choose an 
appropriate self-similar coordinate. 

The energy spectrum in the ejected matter is deduced 
from the solution in terms of the mass M(> e) of the 
ejected matter with the Lorentz factor greater than a 
certain value. The power law index of this mass with re¬ 
spect to A greater than 1/2 indicates that the enhance¬ 
ment of the initial explosion energy will accelerate the 
matter near the leading edge more effectively than that 
in the inner part of the ejected matter because this re¬ 
lation holds only in the outermost layers of the ejecta. 
Otherwise, the total energy would not be conserved. It 
is also expected from this power law index that if the 
explosion energy is not distributed uniformly in the lat¬ 
eral direction, the non-uniformity will be enhanced near 
the leading edge of the ejected matter and result in jets. 
This solution will be useful to test relativistic hydrody- 
namical codes whether they are appropriate to treat the 
passage of a shock wave across the surface of the matter 
to a vacuum. 


This work has been partially supported by the grant 
in aid (16540213) of the Ministry of Education, Science, 
Culture, and Sports in Japan. 


APPENDIX 

ANALYTICAL SOLUTIONS BEFORE SHOCK EMERGENCE 

Differential equations can be analytically integrated after some manipulations. Using a new 

variable defined as ^ = —(/x, equations P|l -(pn |l are transformed into differential equations with respect 

din/ {m — 6)^ + 8m — 4:6 

d/ -I-(Sd-I-m-I- 8 )^-I-4(m-I-1) ’ 

din 5 m^ + 7m — 8S 

d/ _|_ (3,5 ^ _l_ _l_ _l_ X) ’ 

dlnd (m — d)/^ -I- (lOm — 8 d)^ -|- 18 to — lOd 
d/ (^-I-2){/2 + (3(5 + m-I- 8 )^-I-4(to-I-1)}' 

These equations can be integrated to yield 


/ = 


^ - CD 


(i-m)(l + Ci) 


9 = 


4 + ai 

^-Oil 


^ — 02 




1 -I- Cki 


rt(l + C2) 


1 -I- a;2 


^ - 02 




h = 


where constants have been introduced and are defined as 


/ - Oi 

C3(5—m) 

^ - 02 

1 -1- Oi 


1 - 1-02 


1 -I- 02 

C‘4{<5—m 




Ol ,2 = 


— (3d -I- m -I- 8 ) ± (3d -I- m -I- 8 )^ — 16(m -I- 1) 


Cl = 


— (3d -I- TO -I- 8 ) -I- 8(2to — 1)/(to — d) 


Oi — O 2 


independent 
to /, 

(Al) 

(A2) 

(A3) 

(A4) 

(A5) 

(A 6 ) 

(A7) 

(AS) 
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C 2 


—36 — m + 6(1 — 6/m) 

ai — a2 


(A9) 


(ai — 2)[m? 3(5^ — 4(5m — 25 to -I-19(5) -I- 2(m -|- 2){rn? — 36'^ + 26m — 2m) 

2(m -I- 2)(m — (5 )(q!i — 0 ( 2 ) 

—m? + 36"^ — 26m + 2m a 2 — 2 

^ (ai — 2)(m — (5) ^Ofi —2’ 


(AlO) 

(All) 


C 5 = 1 - C 3 - C 4 . (A12) 

Here we have used the initial conditions m. From (EUl-lEll), we can obtain the behaviors of f, g and h in the limit 
of ^ ^ + 00 , i.e., X where the constant terms Ofi and a 2 in these equations are negligible. Thus we obtain 


•f ^5—m „ /•—m u tS—m 

/CX4 ,/^0C4 , 

which are rewritten in terms of x Ets 

/ oc {-x)^, 9 cx (-x)"^, h oc (-x)^- 

These results justify our anticipation for the limiting behavior of the solution (EHl-iEU- 
For the eigenvalue m = (2 — Vs) 6 obtained in ini equations lldTi-Kali are reduced to 

din/ _ 2(2-73)5 

d/ ^ -f 4 -|- 273 (1 -f 6) 

d In 5 73 (73 — 2)5 

d/ ^ -f 4 -|- 273 (1 -l- 5) 

dlnd (2-7^5 [2 -p 73 73(1-f 25) | 

~W “1 + 73(1 + 5) I C + 2 ^^ + 4 + 273(1 + 5)/■ 


(A13) 

(A14) 

(A15) 

(A16) 

(A17) 


These can be integrated to the analytical expression for the numerical solution presented in im By directly substituting 
the eigenvalue of m into equations (IA4ll - (IA6ll . we also obtain the analytical expressions, 


In f = 


( 9 A l>n 

^ + 4 + 2y/3(l + (5) 


3 + 273(1 + 5) 


/ 9a/9 1 XIti 

+ 4 + 2'\/3(l S) 


3 + 273(1 + 5) 


(A18) 


(A19) 


\nh = 


1 + a/S (1 + (5) 


In 1/ + 2| + (273 - 3 ) (1 + 25) In 


+ 4 + 2'\/3(l H“ S) 


3+ 2^3(1+ 5) 


(A20) 


REFERENCES 


Blandford, R.D. McKee, C. F. 1976, Phys. Fluids, 19, 1130 
Colgate, S. A., McKee, C. R., & Blevins, B. 1972, ApJ, 173, L87 
Costa, E., et al. 1997, Nature, 387, 783 

Fields, B. D., Daigne, F., Casse, M., & Vangioni-Flam, E. 2002, 
ApJ, 581, 389 

Galama, T. J., et al. 1998, Nature, 395, 670 

GandePman, G. M. & Frank-Kamenetskii, D. A. 1956, Soviet Phys. 
Dokk, 1,223 

Iwamoto, K., et al. 1998, Nature, 395, 672 
Johnson, M. H. Sz McKee, C. F. 1971, Phys. Rev. D, 3, 858 
Kazhdan, I. M., & Murzina, M. 1992, ApJ, 400, 192 
Kulkarni, S. R. et al. 1998, Nature, 395, 663 

Landau, L. D. Sz Lifshitz, E. M. 1987, Fluid Mechanics, 2nd ed. 
Pergamon 


Marti,J.M. Sz Muller,E. 1994, J.Fluid Mech., 258, 317 

Matzner,C.D. &; McKee,C.F. 1999, ApJ, 510, 379 

Meszaros, P. Sz Rees, M. J. 1997, ApJ, 476, 232 

Nakamura, K. Sz Shigeyama, T. 2004, ApJ, 610, 888 

Perna,R. & Vietri, M. 2002, ApJ, 569, L47 

Rees, M. J. Sz Meszaros, P. 1992, MNRAS, 258, 41P 

Rhoads, J. E. 1997, ApJ, 487, LI 

Sakurai, A. 1960, Commun.Pure Appl.Math.,13,353 

Tan, J. C., Matzner, C. D. Sz McKee, C. F. 2001,ApJ, 551, 946 

van Paradijs, J., et al. 1997, Nature, 386, 686 
























